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Abstract 

A mechanism is presented, in which the mixing of right handed 
heavy Majorana neutrinos creates a CP— asymmetric universe. When 
these Majorana neutrinos subsequently decay more leptons than anti- 
leptons are produced. The lepton asymmetry created by this new 
mechanism can exceed by a few orders of magnitude any lepton asym- 
metry originating from direct decays. The asymmetry is finally con- 
verted into a baryon asymmetry during the electroweak phase transi- 
tion. 
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The generation of the baryon asymmetry in the universe has been dis- 
cussed in many articles (I], |2[ . Two prominant scenarios are the Grand Unified 
Theories |3], |2|] and the production of an asymmetry through extended solu- 
tions of field theories (sphalerons) 0, ||. At this time both scenarios generate 
asymmetries which are small. This comes about because in the former case 
CP— violation is produced by higher order effects and in the latter the tunnel- 
ing rate through potential barriers is small (in addition, if the higgs particles 
are heavier than 80 GeV, then the baryon asymmetry thus generated will be 
completely erased) ||. 

A third scenario includes heavy Majorana neutrinos whose decays gener- 
ate a lepton asymmetry, which later on is converted into the baryon asymme- 
try. By their very nature Majorana neutrinos posses AL = 2 transitions and 
in addition they may have couplings which allow them to decay into the stan- 
dard higgs and leptons, i.e., N — > <f)U~~ ||. In these models the CP— violation 
is introduced through the interference of tree-level with one-loop diagrams 
in the decays of heavy neutrinos, which we shall call e'— type effects (direct 
CP— violation) |], 0. An additional contribution from self energies was in- 
cluded in ref [8-10]. A new aspect was pointed out when it was realised 
that the heavy physical neutrino states are not CP— or lepton-number- 
eigenstates. Therefore as soon as the physical states are formed there is im- 
printed on them a CP— asymmetry and a lepton-asymmetry. This we shall 
call e— type effects (indirect CP— violation). It is a property which appears 
in the eigenstates of the Hamiltonian at an early epoch of thermodynamic 
equilibrium, even before the temperature of the universe falls down to the 
masses of the heavy Majorana particles. 

We demonstrate the phenomenon with a Gedanken experiment. Consider 
a universe consisting of a large p — p collider which produces s— and s— pairs. 
The s and s quarks hadronize into K° and K° mesons whose superpositions 
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are the physical states 



K L ^^x [{l + e)K°±{l-e)W\. 



The probability of finding a \K° > is proportional to |1 + e| 2 and the proba- 
bility for a \K° > is proportional to |1 — e| 2 which are not equal. When the 
particles decay 



the above asymmetry survives as an asymmetry of the detected e + 's and 
e~'s. The electron-positron asymmetry generated is the same for the Kl and 



In this article we point out that a similar situation arises in the for- 
mation and decays of Majorana neutrinos. We adopt the wave function 
formalism to calculate the eigenstates of the Hamiltonian and their decay 
rates. In this formalism we can extend the region of our calculation to the 
case of very small mass differences using degenerate perturbation theory. We 
find that for small mass differences between the two generations the indirect 
CP— violation (e— type) produces a lepton asymmetry much larger than that 
of the e'— type. 

We work in an extension of the standard model where we include one 
heavy right handed Majorana field per generation of light leptons (Ni, i = 
1, 2, 3). These new fields are singlets with respect to the standard model. The 
lagrangian now contains a Majorana mass term and the Yukawa interactions 
of these fields with the light leptons, 



K L — > 7r e T v and K s — > tt e T v 



K s Q- 



int 



E Mi[(N m yN m + N Ri (N Ri ) c ] 



i 



+ E Ki Nm ^£ La + Yl h <*i ha(f> N Ri 



(1) 
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where <p T = (— is the higgs doublet of the standard model, which 
breaks the electroweak symmetry and gives mass to the fermions; l^a are the 
light leptons, h a i are the complex Yukawa couplings and a is the generation 
index. We have adopted the usual convention for charge conjugation : N c = 
CN T '. Without loss of generality we work in a basis in which the Majorana 
mass matrix is real and diagonal with eigenvalues M«. 

The states |iVj > decay only into leptons, while the states \N° > de- 
cay only into antileptons (figure 1). For this reason the states \Ni > and 




Figure 1: iVj and Nf decaying into l La and {IluY 

\N? > have definite lepton numbers and are the appropriate states to de- 
scribe CP— violation in the leptonic sector. They are analogous to the K° 
and K° states. The idea is as follows: Through the presence of the Yukawa 
interactions we obtain one loop corrections to the mass matrix (figure 2), 
such that the corresponding mass eigenfunctions are no longer the |JVj > 
and \N£ > states, but a mixture of them. It is these physical eigenstates 
which evolve in time with a definite frequency. If they are shown to be asym- 
metric linear combinations of the |iVj > and \N£ >'s then we have created 
a CP— asymmetric universe. By asymmetric linear combinations we mean 
that the |iVj > and \N? >'s enter with different complex phases into the 
decomposition of the eigenfunctions. The subsequent decay of these fields 
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will produce the desired lepton asymmetry. 

As a result even if we start with equal numbers of |JVj > and \N? >, they 
will evolve according to the asymmetric eigenstates of definite time develope- 
ment. Since |iVj > and \N? > carry different lepton numbers given by their 
interactions, this means that a lepton asymmetry is established through mix- 
ing before the fields actually decay. Herein lies the main difference between 
our model and the literature. 

For the sake of simplicity we consider two generations of Majorana neu- 
trinos, where the indices % and j take the values 1 and 2. We assume the 
hierarchy M 2 > M x . In the basis (|A^f > |iV 2 c > \N X > \N 2 >) the effective 
Hamiltonian of this model can be written as, 



Once we include the one loop diagram of figure 2, there is an additional 
contribution to the effective Hamiltonian, which introduces CP— violation in 
the mass matrix. We treat the one loop contributions as a small perturbation 
to the tree level Hamiltonian. 



7?°> 



/ Mi \ 

M 2 

M x 

V M 2 / 



(2) 



/ 








(3) 



ttW ff(l) 

ffW rr(l) 
n 12 n 22 



with 





a a 



f dis 
aij 



(4) 





M^h m h* aj + M^h* m h a3 (g: 



a a 



t dis 
aij 



\9%) (5) 
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Figure 2: One loop contributions to the mass matrix 
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and 



(i) _ nr(i) 



W 1 ' = H 
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2M l J2h ai h* m 



(flffi - ( 6 ) 



as can be easily read off from figure 2. The dispersive part g^fj can be 
absorbed in the wave function and mass renormalization. The absorbtive 
part <7^ of the loop integrals is given by, 



< = lis < 7 > 



neglecting terms of order O (m^/p 2 ), O (rn 2 /p 2 ^j with p 2 > M 2 . 

Employing ordinary first order perturbation theory the eigenfunctions of 
the effective Hamiltonian H = W ^ + are found to be 



*1 


> = 




> +1^ > + ai |iV 2 c >) 




> = 


-i(|iVi > +a 2 \N 2 
yj\ 


> — |iVf > -«i|iV 2 c >) 




> = 


-id ^2 > -«i|iVi 
VA/ 


> + |iV 2 c > -a 2 \N{ >) 




> = 


-r=.{\N 2 > -«i|iVi 


> -\N 2 > c +a 2 \N{ >) 



(8) 



with, 



M ± H^ + M 2 HW 12 _ M 1 H( 1 ) 12 + MaifJ? 

ftl ~~ ^72 772 ' 0:2 ~~ 



Mi-Mi 1 Mi-Mi 

where the normalization factor is A/" = 2 + |«i| 2 + |«2| 2 - 

The states > and \^[ > are eigenstates with mass eigenvalues ±(Mi + 
i?n^). They are related by a chiral 75— transformation and correspond to the 
same physical state. The same holds for the states \^ 2 > and \ty' 2 > with 
eigenvalues ±(M 2 + H^). For the rest of our calculation we shall only 
consider l^i > and >■ 
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The asymmetry parameter can be defined as, 



(9) 



which is a measure of the lepton asymmetry generated when the physical 
states j^j > finally decay This can be calculated using, 



T^^i oc Y \hai + a> 2 h a2 \ 2 

a 

= Y [\ h »i\ 2 + |«2| 2 |^2| 2 + 2Re(a 2 h* al h a2 ] 



Vy^ic OC Y\h* a l + a lh 



* |2 
a2\ 



Y [l^il 2 + \®i\ 2 \h a2 \ 2 + 2Re(a 1 h al h* a2 ] 



Y [\ h ^\ 2 + \ai\ 2 \h a i\ 2 ~ 2Re(a l h al h* a2 ) 



* |2 



Y [\ h ^\ 



|«2| 2 |^. 



al I 



2Re(a 2 h* al h 



a2) 



(10) 



In addition to the CP— violating contribution due to the mixing of the 
states \Ni > and \N? >, which we call 5, there is another contribution e' 
coming from the direct CP— violation through the decays of | iVj > and \N£ >, 



> = i(i-^^i:iM 2 M, 



(id 



which has been discussed in the literature extensively || 0]. 

Then it is straightforward to show that the asymmetry parameter consists 
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of the following two parts, 



A = d + 5 



(12) 



The new indirect CP— violation 5, which enters through the mass matrix, is 
given by, 



Re 



2719 c m[^m! 



1 



+ 



1 



^2a I hal 1 2 So I h 



a2\ 



(13) 



where, 



c = --imE(WEW] 



+ 



7T 



^2a |^al| 2 S« 1^ 



<i2 



(14) 



From this expression it is clear that this contribution becomes significant 
when the two mass eigenvalues are close to each other (figure 3). On the 
other hand the perturbation theory used for this expression is valid only for 
\M 1 -M 2 \ > \H\f\ or \H§\ 

To find out the value of 5 in the vicinity of Mi = M 2 we now write 
Mi = M and M 2 = M+r]M and consider the case Mr] < \Hij \ or equivalently 
V < I J2 a h a ih* a j\- Then the total Hamiltonian, 



H = 



( 

M + H u H 12 

V H 12 (M + rjM) + H 22 



M + H u 
H\ 2 





#12 

(M + r]M) + H 22 





(15) 



will have the eigenvalues, 

A 2 « M 2 (l + r ] ) + M{H 11 + H 22 ±c} 



or 



A 



± 



(M + ± V M) + l - {H u + # 22 ± c} 



(16) 
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r — > 



Figure 3: Comparison of the two CP— violating contributions to the lepton 
asymmetry in units of C, where r = Mf/Mf 
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with c = \][Mr] — (Hu — H22)] + (H12 + Hi 2 ) 2 } neglecting terms of order 

* I 

aj I 



M 2 \ J2 a haih^A 3 - The eigenvectors are now given by. 



|^i > = 1 fliVi > +a\N 2 > + \ N i > +®\ N 2 >) 



2(1 + a 



2\ 



|* 2 > = ; 1 (\N 2 > +(3\Nt > +|iV 2 c > +p\N{ >) (17) 
v/2(l + |/3| 2 ) 



with 



H 12 + #12 . #11 -H 22 -riM-c 

a = , p = — , 

H n -H 2 2-vM + c H 12 + H l2 

Note that the states \ipi > and \tp 2 > are not CP— symmetric, because 
under CP— transformations a\Ni > transforms into ot*\N£ >. With this we 
can now calculate the lepton asymmetry for the degenerate scenario defining 
the asymmetry parameter as before. It is now given by, 

,n (181 

if + Gr 6 ) 2 Re 2 (E>*iM 

This expression vanishes for rj —>■ as expected, because there is no detectable 
mixing between two identical particles (figure 4). The matching of the two 
solutions occurs at r] m Da^ai^c*2/5 and thus we obtained an asymmetry 
for small and large mass differences. This is shown in figures 3 and 4. For 
V ~ I J2 a h a ih* a - 1 a big enhancement of the asymmetry parameter 5, by several 
orders of magnitude, is achieved. The enhancement factor is roughly 1/rj. 
For large values of 77 the two distinct contributions e' and 5 become of the 
same order. 

The dynamical evolution of the lepton number in the universe is governed 
by the Boltzmann equations. This is the appropriate tool to describe any 
deviation from thermal equilibrium. The framework and the derivation of 
the Boltzmann equation is reviewed in ||; we adopt the same notation. 
We make the approximations of kinetic equilibrium and Maxwell-Boltzmann 
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Figure 4: Lepton asymmetry in units of C generated for a small mass differ- 
ence, where P = E« h* al h al and r = Mf /Mf . 
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statistics. The out-of-equilibrium conditions occur, when the temperature T 
drops below the mass scale = Mi + H$ , where the inverse decay is 
effectively frozen out. 

The density of the lepton number asymmetry — n\ — n\ c for the left- 
handed leptons can be shown to evolve in time according to the equation: 



f- + 3Hn L = (e' + S 1 )rf i [n^ - n£] - — n%r» - 2n,n L (a\v\) (19) 



where 5\ is the contribution proportional to 1/ J2 a \ h a i\ 2 i n the expression for 
5. The second term on the left side comes from the expansion of the universe, 
where H is the Hubble constant. is the thermally averaged decay rate of 
the l^i) state, n 7 is the usual photon density and the term (cr|v|) describes 
the thermally averaged cross-section of I + <fy < — > l c + <f) scattering. We note 
that the first term of the right side of this equation describes the creation 
of lepton number and is proportional to (e' + S\), while the last two terms 
are responsible for any depletion of lepton number and are coming from the 
inverse decay and the scattering respectively. The density of the ip\ state 
satisfies the Boltzmann equation, 



In order to find a solution to this set of coupled differential equations it 
turns out to be useful to transform to new variables. We introduce the 
dimensionless variable x = M^/T, a particle density per entropy density 
Yi = rii/s and make use of the relation t = x 2 /2H(x = 1). 

In addition we define the parameter K = F^x = l)/H(x = 1) which is 
a measure of the deviation from equilibrium. For K <^ 1 at T we 
are far from thermal equilibrium so that both inverse decays and 2^2 CP 
non-conserving scattering processes are not important and can be safely ig- 
nored. With these simplifications and the above redefinitions the Boltzmann 




dn^ 
dt 



+ 3ifn^ = -r£(rty 1 




(20) 
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equations effectively read : 

=f = -(Y^-YftKa? (21) 

For very large times the solution for Yl has an asymptotic value which is 
approximately given by 

Y L = ^ = -(e' + 5 1 ) (22) 

s g* 

where g* denotes the total number of effectively massless degrees of freedom 
and is of the order of O(10 2 ) for all usual extensions of the standard model. 
The lepton asymmetry thus generated will then be converted to the baryon 



asymmetry of the universe during the electroweak phase transition |12| and 
is approximately given by ~ \ n L- So we have demonstrated that the 
baryon asymmetry is proportional to (e' + Si). This new contribution was 
shown to be at least of the same order as the e' (as shown in figure 3) and 
for small values of r, 5 exceeds e' by several orders of magnitude (figure 4). 
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